Meissner effect without superconductivity from a chiral d-density wave 
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We demonstrate that the formation of a chiral d-density wave (CDDW) state generates a Topolo- 
gical Meissner effect (TME) in the absence of any kind of superconductivity. The TME is identical 
to the usual superconducting Meissner effect but it appears only for magnetic fields perpendicular 
to the plane while it is absent for in plane fields. The observed enhanced diamagnetic signals in 
the non-superconducting pseudogap regime of the cuprates may find an alternative interpretation 
in terms of a TME, originating from a chiral d-density wave pseudogap. 

PACS numbers: 75.20.-g, 71.27.+a, 74.72.-h 



The Meissner effect is considered to be the most direct 
signature of superconductivity [l[ . However, the surpris- 
ing observations of such enhanced diamagnetic signals 
well above the superconducting transition tempera- 
ture in the pseudogap regime of the cuprates 0], con- 
stitute a fascinating puzzle. There are two proposals 
for the nature of this regime that appear to dominate. 
The first, associates the pseudogap with a d x 2_ y 2 density 
wave (DDW) 0,||, also called orbital antiferromagnet 
[E j which normally competes with superconductivity. 
The second associates the pseudogap with spontaneous 
vortex-antivortex unbinding leading to incoherent super- 
conductivity 0] that should persist well above the su- 
perconducting T c . This theory is reminiscent of the well 
known Kosterlitz-Thouless transition 0] . 

The available ARPES [Ic} and STM [U experiments 
cannot differentiate a SC from a density wave (DW) gap, 
and therefore appear somehow incapable in settling di- 
rectly the issue. On the other hand, the unusual Nernst 
effect and most importantly, the enhanced diamagnetic 
signal that accompanies it for a very large temperature 
region above the SC critical temperature 0|, has been 
considered as a major argument in favor of the incoher- 
ent SC scenario. In fact, the enhanced diamagnetism is 
viewed as a signature of the usual Meissner effect associ- 
ated solely with the SC state, and would therefore contra- 
dict the d x i_ y i density wave scenario since no Meissner 
effect was expected in that case Q. 

In this letter we put forward the Topological Meissner 
effect (TME), that results from a chiral d xy + id x 2_ y 2 
density wave (CDDW) state. In fact, the Nernst region 
of the pseudogap regime may well be associated with a 
CDDW. The most intriguing property of a CDDW is 
that parity (V) and time-reversal (T) violation induces 
Chcrn-Simons terms in the effective action of the elec- 
tromagnetic field, providing the possibility of the TME 
and the Spontaneous Quantum Hall effect (SQHE) ear- 
lier discussed [H, EE El EE EE IS El- As we shall 
demonstrate, the TME is described by the same equa- 
tion we find in the usual Meissner effect of a supercon- 



ductor. Though, its origin is radically different. In our 
system we encounter the realization of Parity Anomaly 
EE EE]> with the emerging Chern-Simons terms provi- 
ding a topological mass to the electromagnetic field, in 
a gauge invariant manner [2E [2lJ. Moreover, the pos- 
session of chirality perpendicular to the plane, implies 
that the TME is strongly anisotropic. Particularly, it 
takes place for magnetic fields perpendicular to the plane 
while it is absent for in plane fields, in accordance with 
the experimental observations '2\ . Note finally that a chi- 
ral d-density wave state has also been shown recently [22j 
to explain the experimental results concerning the Polar 
Kerr effect in YBCO (H- 

In order to demonstrate how the TME arises, we shall 
consider the following BCS hamiltonian for the CDDW 
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which describes a d xy +id x 2_ y 2 state characterized by the 
wave- vector Q = (7r,7r), which is commensurate to the 
lattice (fe + 2Q = k). Since spin degrees of freedom do 
not get involved we have considered spinless electrons, so 
that all our results will refer to one spin component. Fur- 
thermore, we use g„ v — (1,-1,-1), k l = k = (k x ,k y ), 

k» = k= q» = q= (q Q: q), fx = 0,1,2, i = 1,2, 

e > 0, h = 1 and we assume that repeated indices are 
summed. In the derivation of the Chern-Simons terms 
we shall restrict ourselves to the zero temperature case 
while necessary extensions to finite temperatures will be 
afterwards performed. In addition, the summation in 
k— space is all over the whole I s * Brillouin zone rather 
than the reduced Brillouin zone. This implies that the 
operators c k and c fc _|_Q do not describe independent de- 
grees of freedom. 

In Eq.||Tl) we have introduced the CDDW order param- 
eter A k = ?yA sin k x sin k y +iA (cos k x — cos k y ) , where A 
is the modulus of the id x 2_ y 2 order parameter, i] defines 
the relative magnitude of the two components and also 
determines the direction of the chirality of the state. The 
chiral character of the state implies the existence of an 
intrinsic angular momentum in fe— space, perpendicular 
to the plane, originating from V — T violation. Specifi- 
cally, the d x 2_ y 2 component violates T as it is imaginary, 



2 



while the d xy component is odd under V in two dimen- 
sions, which is defined as (k x , k y ) — ► (k x , —k y ). 

In order to obtain the total electronic Hamiltonian H, 
we have to add the corresponding kinetic part TL kin ■ For 
the kinetic part we keep only the nearest neighbors hop- 
ping term e fc = —t (coak x + cosk y ) satisfying the nest- 
ing condition e fe _|_Q = — e fc , while we also set the chem- 
ical potential equal to zero. Our approximation can be 
justified by considering that our system is close to half- 
filling. Under these conditions the excitation spectrum 
consists of two bands which are fully gapped leading 
to the topological quantization of the Hall conductance 
[II El El El EE 53, El , which is the coefficient of the 
Chcrn-Simons terms. Omitting the next nearest neigh- 
bors hopping term 5 k = t' cos k x cos k y does not alter 
qualitatively the occurrence of the TME. However, its 
inclusion would destroy the quantization of the Hall con- 
ductance, as in this case, the system is not fully gapped. 
Similar effects would arise in the presence of disorder or 
by including the z-axis hopping term. 

Under this conditions, the total Hamiltonian of the 



J 

k+Q^k+Q 



system becomes TL = \Yl<k e k ( C I 
^A fc c^c fc _|_q + h.c^j . We obtain a compact representa- 
tion of TL by introducing the spinor = -^(c^ c \-^q)^ 
the isospin Pauli matrices r and the vector g k = 

(KeA k , -XmA k , e fe ). This yields TL = £ fe * { g k ■ r * fc . 
The latter indicates that the ground state of the sys- 
tem depends on the orientation of the g vector in isospin 
space. As a result, this hamiltonian supports skyrmion 
solutions which imply the presence of a Chern-Simons 
action (see e.g. (3). 

To reveal the emerging Chern-Simons terms, we 
have to take into account the fluctuations of the 
U(l) gauge field A^. We add to the Hamiltonian 



the term TL„ 
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Sfc ^k+q^k+q,k^fi( 



(2tt) 2 <^k * k+q^ k+q,k^fJ.yi} * k 

I^Ek *l+ q ^^(-9)A(9)* fc! which describes the 
interaction of the gauge field with the electrons. We 
have introduced the paramagnetic interaction vertex 
T k+ q ,k = ~( e i e §k? 9 k ■ T )i where M = 0,1,2 and i = 
1,2. At one-loop level, the effective action S em is given 
by the relation S em = \ f ^ A»{-q)^ v (q)A» {q), 
with the Polarization tensor n^, defined as \i^ v {q) — 

I Ik Tr (S k n,k +q Gk +q n +q ,k) - ^Pekr Pe ^ the two- 
dimensional electronic density (without including spin), 
Tr denotes trace over isospin indices, Q k is the CDDW 
fcrmionic propagator and we have used the abbreviation 
Ik = / ip Sfc- Computing Ii^ v up to linear order in q, 
yields the Chern-Simons action 



S C s = I ^ ^e llvX A»F", 



(2) 



with Fn V — dn,A v — dvAp. The coefficient of the Chern- 
Simons action is the Hall conductance <r xy . It can be 



shown that it is a topological invariant, reflecting the 
existence of a topologically non trivial, V — T violating 
ground state (see e.g. flS|). Using Eq.j2]) we obtain 
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where we have introduced the winding number of the unit 
vector g k = g k /\g k \, 
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which is equal to 2, because the order parameter com- 
ponents are eigenfunctions of the angular momentum in 
k space with eigenvalue I = 2. 

In the case of a perfect gap, the Hall conductance origi- 
nates only from the chirality N of the lower energy band, 
E k = — \g k \, which is fully occupied. In the same time, 
the upper band, E k = +\g k \, is totally empty while it is 
characterized by opposite chirality. Apparently, if both 
bands were equally occupied then <r xy would be equal to 
zero. In the general case, the two bands, have different 
occupation numbers n_ and n + , yielding a non-quantized 
Hall conductance a xy = §-(Ji- — Deviations from 

nesting, disorder or a chemical potential generally lead 
to such an effect. It is desirable to comprehend, even 
crudely, the effect of these parameters on the Hall con- 
ductance and the TME. 

For this purpose we consider that a finite chemical po- 
tential is added to the system. We shall consider that 
its magnitude is of the order of min\g k \. This minimum 
is realized at the points fco = (i§)=t§)j when -q << 1. 
In this case, we may linearize the spectrum about these 
points so to obtain an approximate analytical solution. 
The two energy bands are described by the dispersions 
E k = -n±^m 2 + (y ■ 5k) 2 , with m = min\g k \ = \g k J, 

Vq the velocity at these points and 6k = k—k . If > m 
and fi < 0, hole-pockets arise in the lower band decreas- 
ing the full occupancy from n_ = 1 to n_ = 1 — n ex , with 
n ex the portion of the empty states. On the other hand, 
if /i > m, electron pockets emerge in the upper band ris- 
ing its occupancy from zero. However, if we take into 
consideration that the two bands have opposite chirality, 
it is evident that in both cases, the effect is the same. 
Consequently, a xy (fj,) — a xy (l — n ex ). The portion of the 
empty states will be determined by the area of the el- 
lipses defined by the four hole-pockets. Straightforward 
calculations yield the simple relation 



(/i 2 -m 2 )/27rfA. 



(5) 



We observe that for small values of \fi\, compared to t 
and A, the effect of doping is negligible. 

We are now in position to obtain the equations of mo- 
tion of the gauge field which will allow us to discuss the 
TME in a Hall bar geometry setup. We consider that the 
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FIG. 1: (Color online) (a) Influence of doping on the Topo- 
logical Meissner effect. The relative change of magnetization 
hardly reaches 1% in the presence of a small chemical po- 
tential (A = 20meV, t = 500meV). (b) The magnetic field 
screening as a function of the position on the Hall bar extend- 
ing from — l x to +l x , for different values of the penetration 
depth A over l x . The magnetic field is totally expelled from 
the sample when l x /X >> 1 exactly as in the superconducting 



Hall bar has dimensions L x — 2l x ,L y » L x extending 
from — l x to l x on the x-axis. The relation L y >> L x in- 
dicates that there is negligible y-dependence of the gauge 
fields (d y = 0). To describe the dynamics of the prop- 
agating gauge field we have to add in Eq.© the three- 
dimensional F 2 kinetic term multiplied with the 2-axis 
thickness d. The final gauge field action is 



S em = I d*x <j --F^F^ + -e MvX A*F 



(6) 



where 9(x) — a xy in the bulk of the CDDW, which is 
considered homogeneous. Variation of S em yields 



9"^ + ^^ = -^Va(^M A . (7) 



iX 



(8) 



where n is the unit vector normal to the boundary surface 
a. Eq.||7|) describes the dynamics of while Eq.© 
provides the boundary conditions. In both equations the 
terms on the right hand side stem from the fact that we 
are dealing with a bounded system and the Chern-Simons 
action is not gauge invariant on the boundary surfaces. 
In such cases, gauge invariance is recovered by current 
carrying chiral edge modes 24]. In the rest, we assume 
the such modes do exist and extinguish the right hand 
sides of Eq.©,©. 

Using the Coulomb gauge, V • A = 0, and considering 
the static limit, we obtain the following equations 



dE x 
dx 



IT ^ 



and 



1 dB z 



\x dx 



xy TP 

d 



0, (9) 



FIG. 2: (Color online) Magnetization ratio M/Mo versus 
temperature T for different values of the ratio l x /X. For 
T < T onS et the system is in a chiral d-density wave state. 
In a BCS approximation, these curves are identical to the 
superconducting case. 



where we have included the electric permittivity e and 
the magnetic permeability [i. To obtain the TME, we 
apply a magnetic field of magnitude Bq perpendicu- 
lar to the sample . The corresponding boundary con- 
ditions are B z (±l x ) = Bq. The magnetic field satis- 



fies the differential equation 



/ a 2 



B y = 0, with 



^ = [i\/ £ /^ )/ a xy]d, the zero temperature penetration 
depth. This is indeed the equation we find in the case 
of a superconductor. Notice that in our case only the 
z-component of the magnetic field is involved, imply- 
ing that the TME takes place only for magnetic fields 
perpendicular to the plane. Solving the above equa- 
tion using the aforementioned boundary conditions yields 
B z (x) = B Q [cosh(x/X)}/cosh(l x /X). In Fig.|Q}>) we plot 
the magnetic field versus the ratio of l x /X throughout 
the whole sample. For l x >> A we have almost com- 
plete screening of the magnetic field. Integration of B z 
over the whole sample yields the magnetization Ai = 
A^ [tanh (l x /\)]/ (l x /\), where we have introduced the 
zero temperature magnetization Mo = BqL x L Vi corre- 
sponding to full penetration. One can readily obtain the 
dependence of the magnetization on the chemical poten- 
tial in the l x >> A regime. Using Eq.©, we find — 
Mq\(jj,)/I x = M/{l-n ex ) and SM/M = n ex /{l-n ex ). 
As shown in Fig.(QJi), the relative change of the magne- 
tization due to doping is unimportant for low doping, 
where it stands 5M/M ~ n ex ~ /i 2 /iA. 

We may also obtain the temperature dependence of the 
magnetization by assuming that the penetration depth 
has the usual BCS temperature dependence [A/A(T)] 2 = 
1 — {T/T onset ) . Under this consideration, we obtain the 
magnetization curves shown in Fig. {5} which are identical 
to the ones encountered in the superconducting case as 
expected within this BCS treatment. If we compare the 
magnetization curves of Fig. © with the experimental re- 
sults in the cuprates, we observe that our BCS approxi- 
mation does not provide a fully satisfactory fit. Nonethe- 
less, a strict quantitative comparison calls for an imple- 
mentation of our picture, that would be more adapted to 
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FIG. 3: (Color online) Spontaneous Quantum Hall Effect as 
a signature of the chiral d-density wave. Symmetrical sur- 
face currents running along the y-axis create a magnetic field 
distribution B(±l x ) = ±Bo leading to the spontaneous gen- 
eration of a Hall voltage Vh- 



the cuprate materials. 

A direct experimental verification of the presence of a 
CDDW could be provided by the SQHE (Fig.©). If a 
CDDW is present, a Hall voltage can be generated by 
the sole application of a magnetic field. Specifically, we 
solve Eq.© with boundary conditions E x (±l x ) = and 
B(±l x ) = ±i?o, as in Ref . (lo| . The spontaneously gen- 
erated Hall voltage is, Vh = vB L x [coth (l x /X) — X/l x ] , 
where v = is the velocity of light in the mate- 

rial. If l x » A, V H = vB L x , which implies that the 
applied magnetic field totally transforms to an electric 
field. However, if there is a number of domains of differ- 
ent chiralities in the sample, the SQHE is not an efficient 
probe of the CDDW. Nevertheless, there is an alterna- 
tive route in detecting a CDDW. We refer to the gapless 
chiral edge modes that exist on the boundary surfaces 
separating the bulk from the vacuum. In order to restore 
the gauge invariance of the Chern-Simons terms, these 
current carrying modes appear, constituting a direct in- 
dication of a CDDW 0,LL7|- 



Based on the V — T violation and the possibility of a 
spontaneous electric Hall response via the SQHE we nat- 
urally expect a CDDW to exhibit a spontaneous ther- 
moelectric Hall response (25[, detectable in principle, in 
a Nernst measurement. The unusual Nernst contribu- 
tion has the same origin with the TME and consequently 
they should scale. This implies that the simultaneous 
presence of both enhanced Nernst and diamagnetic sig- 
nals reported in the pseudogap regime is compatible with 
the assumption of a CDDW state. 

In conclusion, we have proposed an alternative way 
of generating a Meissner effect without invoking in any 
manner superconductivity. In our picture, the existence 
of a chiral d-density wave (CDDW) generates the Topo- 
logical Meissner effect (TME) due to V — T violation. 
The direction of the chirality of the CDDW guarantees 
that the TME takes place only for perpendicular to the 
plane magnetic fields, which is in agreement with the 
diamagnetic observations in the pseudogap regime of the 
cuprates. Moreover, a spontaneous thermoelectric re- 
sponse that accompanies the TME is consistent with the 
observed unusual Nernst signal. Note also, that the pres- 
ence of a CDDW is compatible with the recently observed 
quantum oscillations in YBCO (2(| that reported Fermi 
surface pockets in the nodal areas, possibly indicating 
the doubling of the Brillouin zone. As a matter of fact, 
associating a CDDW with the pseudogap regime seems 
quite promising and undoubtedly further theoretical in- 
vestigation should be performed. 
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